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Abstract Expounded are the parameter measurement for narrow resonance and deter- 
mination of corresponding luminosity at e + e~ collider. The detailed theoretical formulas 
are compiled and the crucial experimental effects on observed cross section are taken into 
account. For luminosity determination, the iteration method is put forth which is mainly 
used to separate the interference effect between resonance and non-resonance decays. 
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1 Introduction 

Resonance are a special kind of particles and their study is of great interest and 
importance in the domain of elementary particle physics. Measurements of resonance 
parameters, such as the mass (Mr), total decay width (r t ), partial decay width of fi- 
nal state / (Tf, where / = e,fi,r indicating the e-pair, /x-pair, and r-pair final states, 
respectively), and corresponding branching ratios are fundamental work for high energy 
experimental physics. For resonances of 1 charmonium and bottomnium, such as J /if), 
if)', if;", ■ ■ -, T(nS), (n = 1, 2, 3, 4, ■ • ■) etc, their resonance parameters can be measured by 
scan experiment at e + e~ colliders [l]-[13]. which is one of basic approaches to understand 
resonances in collision experiment. 

For the experiment using scan method, data are taken at several different energy 
points in the vicinity of the resonance to be measured. The minimization technique is 
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usually applied on the estimator which constructed by the difference between the measured 
number of events and the expected number of events. The latter can be obtained by 
theoretical calculation and Monte carlo simulation. Specially, the expected number of 
events for certain final state / at the point with center-of-mass energy W can be obtained 
by the expression 

N f (W,ff) = L(W)-af s (W,ff) , (1) 

where ff is the parameter vector which contains the information of resonance parameters. 
a obs is the experimentally observed cross section (the detailed description refer to section[3]) 
which is the synthetic cross section including resonance part, continuum part, and their 
interference; and also incorporating the effect due to experiment efficiency L is the 
luminosity which can be acquired through several approaches. 
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Figure 1: Differential cross sections for three QED processes: e + e 
and e + e~ — > 77. The center-of-mass energy is 3.5 GeV. 
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In principle any detectable process can be used for luminosity measurement. How- 
ever in order to achieve high precision, one often selects the process which has larger 
cross section and salient characteristic topology experimentally with accurate theoretical 
calculation of the cross section. From these point of view, the QED processes such as 
e + e~ — ► e + e~, e + e~ — > and e + e" — > 77, are most often adopted for luminosity 

measurements [14]. The lowest order of differential cross sections for these processes are 
shown in Fig. [TJ Experimentally, the response of the detector to each of these reactions is 
quite distinct: efficiencies rely on the charged particle tracking (e + e~), calorimetry (e + e~ 
and 77), muon counter and trigger algorithms. The expected theoretical cross 

sections are calculable in quantum electrodynamics; weak interaction effects are negligible 



for charm and B-factories at the level of 1 per mill. [TB] . 



For luminosity measurement, t 
has to be treated with great care^ 



:_he interference effect in the vicinity of resonance peaks 
Such effect not only distorts the cross section in the 
peak region but also shifts the resonance peak position. Especially, when the cross sections 
of resonance and non-resonance processes are compatible, the interference effect are too 
prominent to be neglected. In such circumstance, when we consider how to determine 
the luminosity, we come across a dilemma. On one hand, to determine the luminosity we 
must subtract the contributions due to resonances and corresponding interference effect. 
This can be realized by correct determination of resonance parameters. On the other 
hand, the measurement of resonance parameters depends on the accurate determination 
of luminosity. That is to say the measurement of resonance parameter and determination 
of luminosity are the cause- consequence interdependence. To resolve such an intertwist 
issue, we recourse to an iteration approach which will be expounded in section HI Before 
that, in section [2] and El presented are the formulas for experimentally observed cross 
section which take into account various experimental effects at e + e~ collider such as 
vacuum polarization, initial radiative correction, and beam energy spread. 

2 Cross Section for 1~ resonance 

In this section, we discuss in detail the experimental corrections on cross section and 
provide the analytic expressions for calculation of experimentally observed cross section. 

2.1 Experimental corrections 

The cross section of the resonance process 

e + e- -> Res. -»• / , (2) 
where / denotes a certain kind of final state, is described by the Breit-Wigner formula 

aBw{W) ~ (MjW' (3) 

where W is the center-of-mass energy, T e and Tf are the widths of the resonance decaying 
into e + e~ and /, T and M are the total width and mass of resonance. Taking the initial 
state radiative (ISR) correction into consideration, the cross section becomes [TH] 

<Tr.c(W) = J dxF(x, s) _ * _ a BW (s(l - x)), (4) 





For 77 final state, since only the continuum process exists, there is no interference dilemma and 
luminosity measurement is simple. 



where s = W 2 , x m = l — s'/s, y/s* is the experimentally required minimum invariant mass 
of the final state / after losing energy due to multi-photon emission; F(x, s) has been 
calculated in many references [T6J (TTJ [HI [19] and n(s(l — x)) is the vacuum polarization 
factor. The radiative correction in the final states are usually not considered [201, [21]. The 
reasons are twofold. In the first place, the hadronic final system is very complicated and 
since the radiative corrections depend upon the details of how the experiment is done, it 
is difficult to give a general, model-independent prescription for them. The second reason 
is that our understanding of the hadronic problem is so crude that there is no need to 
worry about the electromagnetic corrections^. 

The e + e~ colliders have finite beam energy spread. The beam energy spread function 
G(W, W) is usually a Gaussian distribution: 

I (W-W 1 ) 2 

G(W, W) = -^^e (5) 

where A is the standard deviation of the Gaussian distribution. It varies with the beam 
energy of the collider. For narrow resonances such as J/ip and ?//, A is usually much 
wider than the resonance intrinsic width. Therefore the beam-spreaded resonance cross 
section is the radiatively corrected Breit-Wigner cross section folded with the energy 
spread function: 

oo 

a G (W) = J dW'a r . c XW')G{W',W), (6) 
o 

where a r , c , is defined by Eq. (@J. 

Take resonance ip' as an example, Fig. [2] displays three cross sections: the Breit-Wigner 
cross section of Eq. Q; the cross section after radiative correction by Eq. and the 
beam-spreaded cross section by Eq. ([6]). From the three curves in Fig. [21 it can be seen 
that the radiative correction reduces the height of the resonance. It also shifts the peak 
position to above the if)(2S) nominal mass. The reduction factor p and the shift of the 
peak A v / i max are approximately expressed by [22] 



where t is defined as 
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2 In any case, if we find later on that it is necessary to do radiative corrections to the hadronic states 
for some specific problem, we can do the calculation then, because the initial state radiative corrections 
and final state radiative corrections can be decoupled to a large extent. 



i 




3.66 



3.67 



3.68 



3.69 



3.7 



3.71 



H„„ (GeV) 



cm 



Figure 2: Inclusive hadronic cross sections for resonance of if)': obw for Breit-Wigner 
cross section; o r c the cross section with radiative correction; o G the beam-spreaded cross 
section with the energy spread effect. In the calculation of these cross sections, the 
following parameters are used [13]: M = 3.68596 GeV, T t = 300 keV, T e = 2.19 keV, and 
A = 1.3 MeV. 
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a ,-k" 1, 

Herein the reduction factor p is the ratio of the maximum resonance cross section with 
radiative correction to that of Born order. 

At the ip' mass, t w 0.0779 and 5 ~ 0.06, then the reduction factor p w 0.51 and the 
shift of the peak &-\fs max ~ 9 keV. The energy spread further lowers down and shifts 
the experimentally observed if)' peak. In the case of a collider with A = 1.3 MeV, the 
maximum height of the ip{2S) peak becomes 640 nb, and the position of the peak is shifted 
by 0.14 MeV above the ip' nominal mass. 

From the example discussed above, the effects due to experimental corrections on 
resonance cross section are fairly prominent. For the continuum, effects are comparatively 
moderate [23]. Next, we will discuss in detail the ISR effect on Breit-Wigner cross section. 



2.2 Initial state radiative correction 

The ISR correction scheme used by earlier experiments, is based on the work of Bon- 
neau and Martin [23] and that of Jackson and Scharre [25J. The former only calculated to 
a 3 order which is insufficient for resonances; while the latter made some mistakes 



The drawbacks due to the treatment of the radiative correction with these two schemes 
were studied for Z in Ref . [26] and for narrow resonances of ip and T families in Ref . [27] . 
In the eighties of the last century Kuraev and Fadin treats ISR correction based on the 
structure function approach which achieves 0.1% accuracy |16j . Afterwards such an ap- 
proach is extensively used for ISR correction which is also followed in this paper. 

The calculation of F(x,s) is summarized in Ref [22]. But for the discussions on 
resonances in this paper, a different form derived in Ref |28j is more useful: 



F(x,s) 
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Here the conversion of soft photons into real e + e~ pairs is included. Notice that x < 1 
and t < 1, so the omitted terms in the above equation are small quantities. Then using 
equality [29] 
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It can be obtained finally 
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In fact, a more simplified formula can be used , viz. 



ay,,, (a) = |(1 + ^"^(cos/M) + (-t + 0(cos/M + 1) j . (17) 

For resonances of ip and T families, the accuracy of expression ( fT3l) is better than 
0.1% while that of expression ( jTTl) is better than 0.2% over a sufficient large energy range 
around the resonances which is usually scanned by the experiments. Therefore even the 
latter is accurate enough to be used for present data fit |28j. 

One remark is in order here. As aforementioned, in Eq. ([3]), T e and Tf are the partial 
widths of the e + e~ mode and the final state / (here / usually indicates the hadronic 
final state) respectively. Here T e describes the coupling strength of the resonance to e + e~ 
through a virtual photon. For example, in potential model, T e is related to the wave 
function at the origin ip(0) in the way 

4a 2 Q>(0)| 2 
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where Q q is the charge carried by the quark in the quarkonium and a is the QED fine 
structure constant. Since the decay of a quarkonium 1 state to e + e~ pair is through a 
virtual photon, there is always vacuum polarization associated with this process. So the 
experimentally observed e + e~ partial width, denoted explicitly as Tl xp , is related to r e 
by the expression 
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This is the convention of Ref. [201 EZ] which is adopted by PDG. In this convention Y e 
means Yf v . So in the above discussion, the factor of vacuum polarization has been ab- 
sorbed into the partial decay width of e + e~ final state, as in Eqs. (|T3l) and ( JTTl) . However, 
when the leptonic decay is concerned, only one polarization factor can be absorbed into 
T e or T^. Therefore in the following formulas, the vacuum polarization factor will be given 
explicitly. 



2.3 Vacuum polarization 

A pedagogical description on the calculation of vacuum polarization can be found in 
many textbooks on quantum electrodynamic, e.g. Ref. [SO]. In this section, we merely 
collected the formulas for the following usage. 

In the actually calculation, the polarization factor is often expressed as 1 + S vac with 
relation 

1 . . 



§ According to the calculation of field theory, 

8 vac = -2-ReU(s) 
-2 • ReU(t) 
-ReU(s) ■ ReU(t) 



(s channel ) 

(t channel ) (19) 
interference bwteen s and t) 



2.3.1 Leptonic part 

When s > 4m 2 , mi is the mass of lepton (I = e, fi, r) 



ReU(s } mf) = - 
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For t < 0, t = — s • sin 2 |(#: polar angle), so Rell(t, mf) is the function of the center-of- 
mass energy E cm and scattering angle 6. 

ReU(t,mf 
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2.3.2 Hadronic part 

For s channel I3T1 
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3 The definition of Il(s) in the literature varries by a minus sign, e.g. between Ref. [16] and 



In the above equation, the summation includes all 1 resonances, such as p, u, 0, J /if), 
t/>', T and so forth. Where R(si) is R- value at In the second term of Eq. ( 1261) . R(si) 
can be expressed through cr/^Si) as 



where &h( s i) indicates the hadronic cross section produced in e + e~ collider at s — si. 
When s < 4ml, °7i( s ) = 0- 

For t channel, the variable of s should be changed into t in Eq. (|26|) . but the value si 
should be kept the same. 

It should be noticed that when a resonance is to be fit, the contribution from itself 
is not included in the summation of Eq. (126|) to calculate of the vacuum polarization. 
In this case, the vacuum polarization is always a smooth function in the vicinity of the 
resonance and can be treated as a Constance. Always the calculation program is readily 
available [32] . 

3 Experimentally observed cross sections 

In e + e~ colliding beam experiment, for a final state /, besides the decays from reso- 
nance which are produced by e + e~ annihilation (refer to Eq. (121)), most often the process 



produces the same final state simultaneously, which is indistinguishable from that due 
to resonance decays. So for the final state produced in e + e~ experiment, it generally 
composes of three parts: the resonance, the continuum, and the interference between 
them. 

In addition, another factor which should be considered is the experimental acceptance 
which actually includes trigger efficiency, reconstruction efficiency, and selection efficiency. 
For the last term, it includes the geometry efficiency implicitly (through the sub-detector 
coverage) or explicitly (by applying certain angle cut). 

Therefore, the so-call experimentally observed cross sections denote total cross sections 
which include interference effect and acceptance. In the content that follows, we will 
present the detailed formulas of cross section for inclusive hadronic final stated, fi + fi~ 
final state and e + e~ final state respectively. 

4 Si is a large energy scale, e.g. in the program of Berends Si — (9.5 2 + 10) GeV 2 . 

5 For exclusive process, form factor has usually to be taken into account, see the details in Refs. [351 
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3.1 Hadronic final state 



The experimentally observed cross section for inclusive hadronic final state is as follows 

af{W) = A«a«(W) + A c h a c h (W) . (29) 

§ Herein the observed cross section has taken into account the effect due to the energy- 
spread, that is 

oo 

af CJ (W) = J dW'G{W,W')af CJ {W) , (30) 
o 

where and A% are acceptances for resonance and continuum hadronic events respec- 
tively [8j; / indicates the final state which can be inclusive hadron (h), /i-pair (//), and 
e-pair (e) respectively; R, C, and I denote resonance, continuum, and interference re- 
spectively. In this section, a denotes the ISR corrected cross section (which is denoted by 
ov.c. i n the previous section). 

In Eq. ( )29|) . the resonance cross section has been given in Eq. ( jTTI) ; for non-resonance 
part: 

o° = Ra°(QED) = = (31) 

h ^ ' 3W 2 W 2 (GeV 2 ) 

where R is the .R-value with radiative correction, viz. 

R = R-{l + 5 rad ) (32) 
where 5 ra d is radiative correction factor at the non-resonance region 



3.2 fi + /i final state 

The experimentally observed cross section for fi + fi~ final state is 

a o bs{w) = A R . a R {w) + A C . a C {w) + A / . ^ {W) (33) 

As indicated in Eq. §6§ the effect due to the energy spread has been taken into account 
for the cross section in cr° bs . The a, which is denoted as <J r , c ,(W) in Eq. ([6]), can be 
approximated analytically by the following form 
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6 For the inclusive hadronic final state, a 1 cx Qii where Q, t is the charge of the quark flavor. So for 
charmonium region, 3a 1 after summation of u, d and s quarks. 
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where s m = 4:E 2 ut , E cut is the energy cut for /i, with E cut > m M . The definition of the 
other variables C\ and C 2 read 
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In the above expression, the terms with yT^T^ indicate the interference part while the 
terms with r e r M are for the resonance part. The variables t, and a are given in Eqs. (jHJ), 
( TTBl and ( |T4|) respectively 



3.3 e + e final state 

The experimentally observed cross section for e + e~ final state is 

a f s (W, ti) = A*($)a*{W) + Af ED (^)af(W) + A{(<d)al nt (W) (40) 

Since QED cross section of e + e~ final state is divergent at the small angle, the acceptance 
of A{§) is relevant to the certain Monte Carlo simulation angl^], that is the e + e~ events 
are produced within the scope ($ — > n — 

The special expressions for the cross section in the above equation are as follows: 
1) for resonance 

0*0*0 = — (i + <W«))- 

s 

7 It should noted that the event produced angle | cos"d prt i\ (p: produce) must be greater than the event 
selection angle |cos$ se ;| (sel: selection). 



i i 



|t ro (1 + 5)a t ~ 2 <P(cosP,t) + [T m 6t + T m t(l + §)]a t-i 0( cos ^+ *) j . (41) 
2) for continuum or QED process 
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3) for the interference 
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The special meaning of parameters in Eq. (I43p is 
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with the variables t and given in Eqs. (Q and (|T6|) respectively. In addition, 
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4 Fit of experimental data 



Assume that the scan data are taken at n pt points with different energies Wi (i = 
1, 2, • • • , n pt ), at each point, the number of experimentally selected events for final state / 
is denoted as n l j (/ = /i(hdaron), e(e-pair), /x(/x-pair)); the theoretically expected number 
is calculated by the formula 

N}(ff) = V ■ a}(ff) , (44) 

where L 1 is the luminosity measured at energy Wi and the energy dependence has been 
denoted simply by the subscribe "i" ; a indicates the observed cross sections given in the 
previous section with the superscription "o6s" removed in this section; ff is the parameter 
vector which contains all parameters to be fit, such as resonance mass (Mr), total decay 
width (r t ), partial decay width of final state / (Tf), energy spread (A), and so on. For 
example, for ip' scan (assuming e-// universality, i.e. T e = T /t ), 

ff=ri{M R ,r t ,r^A) . (45) 

Then chi-square estimator [38] can be constructed as follow^ 

Minimizing x 2 yields best values (estimates) of the parameters wanted [38]. Usually 
Poisson distribution is assumed for the data, then the relation 5n^ = y^} is always 
adopted for the statistical uncertainty of the data. 

As aforementioned, for physics analysis the luminosity is often determined by some 
physics process with salient characteristic topology and large cross section. As an example, 
for ip' scan the luminosity is calculated by e + e" event as follows 

r< _ nf (Wi) (4?) 



However, as mentioned before, among e + e~ events besides the contribution due to QED 
process there is also the contribution from the resonance decay and interference. To know 
n^(Wi), resonance parameters must be determined first. To solve this intertwist difficulty, 
an iteration method is adopted. 

For the j-th iteration, the luminosity at energy Wi is calculated as follows: 



n e {Wj) 

aim 



H» = ^ , (48) 



8 Here for briefness, we only consider the uncorrelated form of estimator. More complicated form with 
the correlation between data taking into account could be referred to Refs. [551 HOI HTj . 



where n e (Wi) is the observed number of events of e + e~ final state, and a{{Wi) is the 
observed cross section (with superscript "ofes" removed) calculated by Eq. (IIP!) at j-th 
iteration. As the first step in the fit, the o\(Wi) is obtained by a rough guess of the 
parameters, or the previously measured ones, if they are available {e.g. PDG values), 
then calculated by Eq. ( 1481) is the Ufa which can be used to work out the expected 
numbers of events for the processes interested, such as N l h , (refer to Eq. (1441)). Then 
utilize the estimator of Eq. (146]) to get the fit parameters 

^ (1) = ^(M« ) r«,rW ) A( 1 )). 

At the next step of the fit, the estimated resonance parameter values can be obtained 
directly from the measured data, which in turn are used to acquire the approximate total 
observed e + e~ cross section <7g(Wj) by Eq. (liOjl . 
With the recalculated luminosity, the parameters 

are fitted again. Such a recursive iteration is repeatly carried out until the corrected 
values of Uu\ are converged in two successive iterations, that is 

Here a is the convergence precision which can be adjusted to meet the need of scan fit. 

5 Summary 

For e + e~ collision experiment, the approximate analytic formulas of cross sections are 
presented for final states of inclusive hadron, e-pair, and /i-pair, where the experimental 
effects are taken into account including initial radiative correction, vacuum polarization, 
and energy spread. The experimentally observed cross sections are also presented which 
take into account of the acceptance of events and the contributions of the resonance, the 
continuum, and the interference between them. 

In the light of e-fi universality, the iteration technique is adopted to figure out the 
cause-consequence interdependence between the measurement of resonance parameter and 
determination of luminosity. Such a kind of methods have been used successfully for 
resonance parameters measurement of ip' by BES collaboration [8]. 
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